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Abstract 

In this paper 1 we consider the Skorokhod embedding problem in Brownian motion. 
In particular, we give a solution based on the local time at zero of a variably skewed 
Brownian motion related to the underlying Brownian motion. Special cases of the 
construction include the Azema-Yor and Vallois embeddings. In turn, the construction 
has an interpretation in the Chacon-Walsh framework. 



1 Introduction 



The S korokhod embedding problem was first proposed, and then solved, by Skorokhodl 
l|l965|) . and may be described thus: 

Given a Brownian motion (B t )t>o an d a centred target law /i can we find a 
stopping time T such that Bt has distribution pP. 

Skorokhod gave an explicit construction of the stopping time T in terms of independent 
random variables, and in doing so showed that any zero- mean probability measure may be 
embedded in Brownia n motion . Since the problem was posed many more solutions have been 
given, see for ex ample 



survey article of |Obl6i 



Dubinsl 1 1968^1 . Root 1 1969^ and Bassl 1 19831) . and the comprehensive 



In this paper we are concerned with t he co nnection betw een two of the well-known con- 
Azema and Yoi 1 1979a ) and Valloisl 1 19831) , and the relationship between 



structions due to 
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these c onstructions and the construction of 
(1979a) construction defines a stopping time 



Chacon and Walsl 



hHHzl) 



. The 



Azema and Yor 



T 



AY 



inf{< > : B t < ip(S t )} 



where St is the maximum of the Brownian motion B and ip is an increasing function depen- 
dent on li. This is an embedding of [i and has a nice optimality property; namely that it 
maximises the distibution of the maxim um for embeddings under which B t ^T is uniformly 
integrable (see Azema and Yor |j^7 9bj)). The other well-known construction of interest to 
us here is the embedding of Ivallois ( lQSflj) . This construction uses the local time Lf of the 
Brownian motion at 0. In particular, for certain functions a, f3, 

T v = w£{t>0:Bt?((3(L?),a(L?))} 

is an embedding. In Valloisl 1 19921) it is shown that for a convex function \1/ a specific choice 



of a and (3 minimises E\&(L|?) among the class of all embeddings of /i. Moreover the same 
choice of a and j3 is optimal for each It is this specific choice of a and /?, for which a is 
decreasing and (3 increasing, which interests us here. 

It is not surprising given Levy's result on the equivalence in law between (St, St — B t ) and 
(Lf , \ B t \) that there is a strong connection between the two constructions (not least when /i 
is symmetric) ; however the aim of this paper is to demonstrate that both embeddings can be 
seen as examples of a much broader class of embeddings. In particular, our embeddings will 
be characterised by a functi on G : R + — > R, and r e lated to variably skewed Brownian motion. 
Skewed Brownian motion ( Harrison and Sheprl 1981 ) can be thought of as a Brownian 
motion with excursions being positive with probability p and negative with probability 



1 — p; for a variably skewed Brownian motion ijBarlow et al.ll200fj|) the probability of positive 



excursions will vary with the local time at zero of the process. 

Given a Brownian motion B t , a variably skewed Brownian motion X t is the solution to 



X t = Bt- G(L t 



(1) 



where Lf is the local time of X t at and G is a function with Lipschitz constant at most 
1 . Our embedding will then be a stopping time of the form 



T = mf{t >Q:X t ? (b(Lf ) , a(L? ))} = inf{t > : B t $ (p(Lf),a(Lf))} 



(2) 



for suitable pairs of functions (a, 6) and (a, 0). When G = 0, X t is just B t (this remark can 
be restated in a more complicated fashion by saying that X is a skew-Brownian motion with 
constant skewness parameter 1/2) and the stopping time is the same as Vallois' stopping 
time. When G(x) = x for all x, we find that X t has only negative excursions and L 
sup s<t B s ; in particular, the embedding is now of the form 

T = inf{t > : B t < */>{&))} 



x 



for some ip and we recover the Azema- Yor stopping time. 
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The embeddings we cons truct have an i nterpretat i on wi thin the balayage construction 



Chacon and Walshl (jl9T6h . As shown in Meiliisor] l|l983|) . the Azema-Yor construction 



of 

is already known to fit in this framework. We shall show how our construction can be 
interpreted in a similar way. 

The paper proceeds in the following way: we shall begin by establishing the important 
quantities in the construction, and some relationships between them; then we show in The- 
orem 13 that there is a strong solution to and hence that given the Brownian motion, 
we may indeed construct X. The main result of this paper is Theorem 1111 where we use 
excursion arguments to confirm that the stopping time is an embedding, and that B t ^r is 
uniformly intcgrable. In Section 0] we discuss the connection with the Azema-Yor and Val- 
lois embeddings, and describe a new embedding in the general framework which maximises 
E , 5(Ly) for 'J convex and the local time of the Brownian motion at a non-zero level x. 

2 The Construction 

Our aim is to show that we can construct embeddings of a target distribution [i in a Brownian 
motion (B t )t>o with Bq = 0, where the embedding is determined by a curve 7 : [0, C] — > 
R x M + , 7(s) = (F(s),h(s)) and what is in some sense the local time of the Brownian 
motion on the function F. The parameterisation of the curve is essentially arbitrary, but for 
simplicity we might choose s to be the arc length . The c onstr uction we prov ide will have as 



special cases the embeddings of lAzema and Yorl l|1979al) and IVallo is (1983). In particular 



we will see that the Azema-Yor construction corresponds to the choice F(s) — h(s) and the 
Vallois construction to the choice F = 0. 

In order to perform the construction, we need to make some assumptions on the behaviour 
of 7. The exact form of the conditions is dependent on the target distribution, and further 
notation is required to state them clearly, however the assumptions have a straight-forward 
graphical interpretation. 

As with many embeddings, there are strong connections with potential theory. We define 
the function 

c(x) = f \y-x\ti{dy)=W\X-x\. (3) 

This is related to a common definition of the potential of a Brownian motio n (in 1-dimcnsion) 
via c{x) = —u^x) where is the potential of /j, (see e.g. IChaco n (1977)). 



Proposition 1. Let fi be an integrable, centred target distribution. Then the function c as 
defined in (J3J) has the following useful properties: 

(i) c is a positive, continuous, convex function such that c(x) > \x\; 

(ii) as x —>■ ±00, c{x) — \x\ — > 0; 

(in) c has left and right derivatives, given by: 

c'_(x) = 2n((-oo,x)) - 1, c' + (x) = 2fj,((-oo,x]) - 1, 
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Figure 1: The 'potential' c(x) for a distribution \x. The embedding is determined by the 
curve 7 = (F, h) and 9(s) and 4*{s) are the x-coordinates of the points at which the tangents 
to c intersect 7 at 7(s). 



so in particular, c is differentiable Lebesgue-almost everywhere. The set {x : fi({x}) > 
0} is precisely the set where c is not differentiable. 

We now discuss the set of admissible curves 7. We suppose that the functions F 
and h are both absolutely continuous, and hence have derivatives Lebesgue-almost every- 
where. We shall assume for convenience that the parameterisation is such that s is the 
arc-length, however this is not necessary, and the extension to an arbitrary parameteri- 
sation is immediate provided that Leb({s : h'(s) — F'(s) — 0}) = 0. We also define 
( = sup{s > : s < c(F(s))}, and require that the curve 7(s) = (F(s), h(s)) hits c, at least 
in the limit s — > £ (although this could occur at ±00). 

The curve 7 is key to defining the other primary variables. The functions 8, (j> are defined 
in terms of c, F and h so that for s > 0, 

6(s) = sup Li: c(A) > Ms) + (A - F(s)) C ^ "1^ VA < x] , (4) 
I x-F(s) J 



©(*) = inf jr e R : c{\) > /»(«) - (A - /•> — VA ■ .< \ . (. r >) 



h(s) — c(x) 

Alternative definitions of 9 and <f> are given by 



9{s) = argminj ^ j . 

Hf)_ 

x<F{s) I F(s) 



(j)(s) = argmax <^ v / \ , (7) 



4 



but note that the argmin and argmax in (jSJ and J7J) may not be uniquely defined. There 
are no such problems with the definitions in Q and J5|. 

In order for the construction to make sense, we must restrict 7 so that for any particular 
value of s the future of the curve will lie inside the (closed) area constructed via the tangents 
at 9(s) and <j>(s) and the potential c. We note that under this restriction, the curve can have 
h'(s) < 0, which is why we cannot paramatrise the curve by its height. Our first aim is to 
express this idea mathematically. Define 

c(9(s))~h(s) . c(x) - h(s) 
R ( s ) = 9(s)-F(s) = x >%) x-F(s) ' (8) 
h(s) - c((j)(s)) h(s) - c(x) 

F(s)-<p(s) x<F (s) F(s)~x 

r W . MzM. m 

In particular, R(s) and S(s) are the gradients of tangents to c which hit 7 = (F, h) at 
(F(s), h(s)). As well as T representing half the difference in the slopes it will also have an 
interpretation as the probability that the embedding has not yet occured. If we use the 
notation f, to denote the right derivative of / then the condition on 7 is 

h' + (s) > max(F! h (s)S(s),F! h (s)R(s)). (11) 

The construction of the stopping time can be considered via the transformations shown 
in Figure El We want to be able to make sense of the 'local time' of B t on F, which we 
denote by Y t . For the present we are not specific about the process Y except that it is an 
increasing process which only grows when is on i 7 " in a sense which will become clear. 
The full definition of Y can be obtained from reversing the construction process we are 
about to describe. 

The curves in the left graph of Figure El (starting with the highest) represent 6, F and <f>. 
In order to define the candidate stopping time we need to transform the horizontal and verti- 
cal scales of this graph. Firstly, we map the horizontal axis by H^ 1 (where H will be defined 
later), so that the curves in the centre graph are plots of F(H~ 1 (l)), 9(H~ 1 (l)), 0(_ff _1 (Z)) = 
(G(l), a(l),f3(l)). Secondly, we transform the vertical scale by subtracting G(l) from the ver- 
tical co-ordinate, so that the third graph shows a(l) — a(l) — G(l) and b(l) = (3(1) — G(l). 
Also shown in Figure El (by the vertical lines) are the excusions of the Brownian motion 
away from F. In the first and second graph these excursions are plotted in the co-ordinates 
(Yt, B t ) and (Lf , B t ) respectively. In the final plot we have a process X t making excursions 
from 0, and it now makes sense to define the horizontal axis to be the (symmetric) local 
time of X t at 0. In particular, if the process in the final graph has local time Lf, then X 
must satisfy 

X t =B t -G(L*). 

We will later show (Theorem EJ) that this equation has a strong solution. The candidate 
embedding is then of the form 

T = M{t>0:X t ?(b(L?),a(L?))} 
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Figure 2: Three different interpretations of the stopping time; the first graph shows the 
functions 9(s),cf>(s) and F(s); the horizontal axis represents the current 'position along the 
curve' or the 'local time' of B on F labelled by Y t — our primary aim is to understand how 
to interpret this. In the second graph we transform the horizontal scale using a function 
H to be determined, so that the horizontal axis is the local time Lf. For the final graph 
we transform the vertical scale by subtracting G — F o H~ x , so that the vertical axis is 
X t — B t — G(Lf). The vertical lines are the excursions of B from F, or excursions of X 
from zero. In each interpretation the stopping time can be though of as the first time an 
excursion goes outside the region bounded by the upper and lower curves. 



where a(l) = fl^" 1 ^)) - F^ff-^)) and 6(0 = ^(H^Q)) - i^if" 1 ^)). This stopping rule 
can then be re-expressed in the form (J2J . 

Having described the form and philosophy of the construction it remains to determine 
the transformation function H, and then to show that the algorithm outlined above leads 
to a genuine embedding of the target law /i. The advantage of representing the construction 
via th e diagrams in Figur es d and is to show how these representations are special cases 
of the Chacon and Walsh] 1 19761) embedding. 

The conditions on the curve 7 ensure that the function 9 is non-increasing, and <j) is 
non-decreasing. It also follows from the definitions that 9 and <j) are both monotonic and 
left-continuous. Note that a and /3 share these properties, but that a and b need not be 
monotonic. 

The first results discuss the properties of the functions R and S introduced in (JSJ) and 
©. Note that (by Proposition ^iii)) we have the inequalities 

l-2/i([0(s),oo)) <R(s)< l-2/i((0(s),oc)), (12) 
2/i((-oo, #«)))-! <S(s)< 2m((-oo, #*)])-!. (13) 
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Lemma 2. The functions R and S are absolutely continuous on closed subsets of [0, £) and 

R( S ) = [^^ S )-n^w]^ 

= £^ y ^ y [/ l '( S )-F'(,)5(,)]^. 

Proof. The function i?(s) is the gradient of the tangent to c that hits 7 at (F(s),h(s)). It 
is clear that i? is non-increasing (by the assumptions on 7), and for every gradient there is 
a unique tangent which hits 7 at some point; hence R must be continuous. 

Also, for 5 small (we insist that F(s — S) < 6(s)) we have F(s — 8) < F(s) + AS and 
h(s - 5) > h(s) - BS for some A, B > 0, and by © 
c(9(s)) - h(s) + B5 



R(s -S) < 



8(s) - F{s) - AS 
c(0(s))-h(s) c(9(s))-h{s) c(6(s))-h(s) BS 



s) - F(s) 9(s) - F(s) - AS 6{s) - F(s) 9(s) - F(s) - AS 

S 



A c(6(s)) - h(s) B 



(s) - F{s) - BS 



e(s)-F( s ) 

so that R is locally Lipschitz on [0, () and absolutely continuous on [0, z] for z < £. In 
particular it is differentiable Lebesgue-almost everywhere, and if R' is a version of the 
derivative of R, then we must have 

R(s) = l+( R'{z)dz 



for s < C (and hence by continuity, for s = £). 

We can calculate the left-derivative of R(s) , using the fact that 6 is left continuous and 
decreasing, to be: 



d~ ( h'_(s) , FL(s)R(s) , 9>_(s) 



ds v ' 6(s) - F(s) 6(s) - F{s) 0(s) - F(s) 



c , Ws)) _ cjmi - h (*) 



However the second term will disappear: if there is no atom at 9(s), it is easy to see that 

while if there is an atom at 0(s), 9'_(s) = (unless s is the largest s corresponding to the 
atom, however then 114|) still holds). 
Consequently 

rfl fa'- 00 FL(s)R(s) 

ds [S> 9(s) - F(s) 9(s) - F{s) 
and the result holds. We can deduce the statement for S(s) similarly. □ 

Corollary 3. The function T is absolutely continuous on closed subsets of [0, Q and Lebesgue- 
almost everywhere differentiable with 

£ = m^m (h ' (s) - F ' m{s)) - mshw) ih ' is) ~ F ' is)s{s)) 

when the derivative exists. 
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Proposition 4. Define 

m _ rWto-rrugytWi, (16) 

Jo r(z) 7 R(z) - S(z) 

[2h'(z) - F'{z) (R(z) + S(z))} R{s) d * s{s y (16) 

Then H is strictly increasing and continuous on [0,C), with continuous inverse H . 
Proof. This follows from and the fact that R(z) = S(z) only if z = £. □ 
Proposition 5. The function G(l) = F o H~ l {l) is absolutely continuous with \G'(l)\ < 1. 
Proof. Both F and H" 1 are differentiable almost everywhere, so that 

[> H>(H-i(l)Y 

It is sufficient to show that |G'(i)| < 1> or equivalently that — JT' < F' < H'. However: 

H' — F' = - ^4 = - F'iJ) > 0, 

R — S R — S R — S 

and / / 

by (dU). □ 

The following lemma follows immediately from Lemma |21 and Proposition 
Lemma 6. We have the relationships 

H' -F' _ R' H' + F' _ S' 

e - f ~ ~ T - ' f-4> ~ T 7 ' 
3 The verification. 

We have now completed the definitions of the auxilliary functions. It remains to show, 
firstly that the process given by is well-defined, and secondly that the stopping time in 
(0) embeds the law \x. 

Theorem 7. Consider the equation 

X t =Bt- G(L?) (17) 

where B t is a Brownian motion, and is the local time at of the process X . Then there 
exists a strong solution to (|17l) . Further, the solution is a variably skewed Brownian motion, 
with probability of positive excursions p given by 

l-G'(l) 



Pd) 



when the local time hf = I. 
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Barlow et al 



2000) under the slightly 



The existence of a strong solution is shown in 
stronger condition that |G'(Z)| < 1. They prove that the solutions of the equation are 
pathwisc unique and that there exists a weak solution, which, by the Yamada-Watanabe 
Theorem, implies the existence of a strong s olution. (The form of the Yamada-Watanabe 
Theorem given in e.g. iRevuz and Yorl 1)199 93 [IX 1.7] is stated for an ordinary SDE, but as 
Barlow et al note, ca n be extended to cover the case of interest.) The proof of Theorem 3.1 in 
carries over w thout alteration to the case where |G'(Z)| < 1, so it only 



Barlow et al 



remains to show the existence of a weak solution in the more general case. The statement 
about the probability of excursions follows on consideration of the case where G(x) = (3x 
for \/3\ < 1. 

Consequently, Theorem will follow from the following result: 
Lemma 8. There exists a weak solution to Q17|l. 



Barlow et al 



The proof of this lemma will follow the proof of Proposition 2.2 of 
the major difference is that it is necessary to use two functions n, in instead of the single 
function <fi used by Barlow et al. 

Proof. For x > define v(x) to be the solution of 

dv _ ( 2 
dx ~ \l + (G'(x))< 

with v(0) = 0, so that 1 < v'(x) < \/2. Define y to be the inverse of v. Then -X < y'(x) < 1. 
Define also 



n(x) 
m(x) 



y{x)-G{y{x)) 
2 

y(x) + G(y(x)) 



(18) 
(19) 



In particular, < n' < y' < 1 and < ml < y' < 1. 
Now define 

\2n'(l)x :x>0 
R(x,l) = { 

\2m'{i)x :x<Q 

Let B t , L t be a Brownian motion and its local time. By applying Revuz and Yor 1 1^99) [VI. 4. 3] 
separately to R+{B u L t ) = 2n'(L t )(B t V 0) and R-{B t ,L t ) = m'(L t )(B t A 0) we deduce 



R(Bt,Lt) - 2 (n'(L s )l {Sam + m'(L s )l {§3<0} ) dB s + ^ (n'(L s ) - m'(L s )) dL s 
= 2/ (n\L s )l {S ^ y+m'(L s )l { ^ <0} )dB s + n(L t )-rn{L t ), 

J 

\R{B u L t )\ = 2^ (n'(L s )l {Ss > 0} - m'(L s )l {§e<0} ) dB s + J (n'(L s ) +m'(L a ))dL s 
= 2/ (n'(L s )l { 3 a > j - m'(L s )l0 a<Q y)dB s + n(L t ) + m(L t ). 
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The stochastic integrals in both of the above equations have quadratic variation 

A t = 4 f (n'(L s ) 2 l {§ m +m / (L s ) 2 l {i3s<0} )d S . 
•> o 

Since n' and m! are bounded, A t < oo a.s., and because n' V m! > ^y' > we also 
have A t | oo a.s.. Consequently, defining T t = inf{s : A s > <}, we deduce that 

Wt = 2 [ \n'(L s )l {Sam +m'(L s )l {Ss<0} )dB s , 
Jo 

) dB„ 



W t = 2| 4 (n'(L s )l fe0} -m'(I s )l { s s<0} ) l 



are Brownian motions. 

Since L t = n(LT t ) + m(LT t ) is increasing only on the set {t : R(BT t ,L,T t ) = 0} and 
|i?(£>T t 7 i>T t )\ — Wt + L t , L t must be the symmetric local time at of R(Br t ,LT t )- 

However, by l(T%|) and l(T§)) . we have 

R(B Tt ,L Tt )=W t -G(L t ), 
and hence a weak solution to (|17(l . □ 

The important fact is that X t is a strong solution, so that it is adapted to the filtration 
Tf of the Brownian motion — in particular, given the Brownian motion, it is possible to 
construct X t and its local time , even if we do not have an explicit construction (see 
however Remark 3.5 of Barlow et al.). 

Given a curve 7, we can define and (j> as in (0J and J5J, and thence define R, S, H and 
G = FoH~ 1 . Define also 

a(l) = 9(H-Hl))-F(H-Hl)) (20) 
6(0 = <l>(H-\l))-F(H-\l)) (21) 

and let X be given by the solution to (|17J) . 

Definition 9. Define the stopping time 

T = T 7 = inf{i > : X t ? (b(L?),a(L?))}. (22) 

T will be our embedding of /i. 

The key idea will be that we can use excursion theory to calculate the distribution of 
the stopped process Bt- We shall first need the following result: 

Proposition 10. Let Y t = iJ _1 (Lf-). The the law ofYr is given by 

¥(Y T > s) =T(s). (23) 
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Proof. By an excursion theory argument, the rate r(l) at which positive excursions from 
zero hit a(l) is the rate at which an excursion from has maximum modulus a(l) multiplied 
by the probability the excursion is positive: 



r(l ) -±_ x 1 - - 
Similarly, the rate of negative excursions is 



l-G'(l) 



2(9(H-i(l))-F(H-i(l))Y 



1 



G'{1) + 1 



G'(Z) + 1 



IKOI 



2(F(H-i(l))-<KH-Hl))Y 



Consequently 
P{Y T > s) = 



T > H(s)) 

H(s) 



exp ■ 



exp 



= exp 

= m 





1 

2 

S T'(z) 

o r(z) 



1-G'(l) 



G'(l) + 1 



2(e(H-Hl))-F(H-i(l))) 2(F(H-i(l))-d>(H-i(l))) 
H'(z)-F'(z) F'{z) + H'{z) 



dl 



9(z) - F(z) 
dz 



F(z) - <j>[z) 



dz 



where the penultimate equality is a consequence of Lemma |H| □ 

Theorem 11. The stopping time T is a UI- embedding: that is Bt ~ /i and the process 
B t AT is uniformly integrable. 

Remar k 12. W e note that T being a Ul-embedding is equivalent to T being minimal in the 
sense of 



Monroe] l)1972j) : that is T is minimal if whenever S is a stopping time with S < T 



and B g ~ u we ha ve S = T a.s.. For further discussions on minimality, see Cox and Hobsoa 



Cox 



Proof. We first show that T is an embedding. Let i be a point such that x > F(Q. By 
definition, 

F(B T >x) = F(X T + G{L%) > x). 

We will only stop at or above x if we have an excursion of X t above a{L^) before the local 
time is too great; in particular, before 6(H^ 1 (l)) < x. When the local time Lf = I, the rate 
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of excursions of X which hit a(l) is p{l)a{l) 1 . Hence, using (J2UJ, (|2U[) and l|23|) . 

?t > a?) = 



1 At _ 

2 ^ WU)) „, 



ny T >zy . 



H\z)dz 



H'(z) - F'(z) 
{z:6(z) > x} 2{9{z)-F{z)) 
_R'{z) 

{z-.e(z)>x} 



T(z) dz 



■ dz. 



Since R is continuous and decreasing, it will attain the lower bound in (1 1 21) on the set 

{z : 9(z) > x}, so that we deduce (since R(0) = 1) 



\B T > x) 



l-(l-2 M ([z,co))) 



A similar calculation shows that for x < 

V(B T < x) = n{{-oo, x]) 

and therefore that Bt ~ M- 

To demonstrate uniform int egrability, we show th at xV(T > H T ) — > as |a;| — > oo, where 



H x = inf{t > : £? t = a;} ("see lAzema et all 1)198(1) or ICoxl (|2£)05)). We consider the case 
x > 0. If either F(s) = h(s) or F(s) = —h(s) then we are reduced to the case of the 



Azema and Yor 1 1979a|) embedding (see Remark [TUl and uniform integrabity follows from 



standard results. Otherwise, the curve 7 enters the interior of the region between \x\ and 
c(x) in FigureHand then sup s< ^ |-F(s)| < K for some K. Then, for x > 2K, 



P(T >H X )= f ¥(L T > I). 

./{(^(h- 1 (;))>#} ■ 



p(0 



■dl. 



x-F{H- x (l)) 

But P(L T > 1) < 1, p(l) < 1 and [x - i^if" 1 ^))]" 1 < (x - K)- 1 < 2/x and so 

xP(T > H x ) < 2Lcb{{Z : 9(H-\l)) > x} = 2H(0- 1 (x)), 
and as x | 00, d^ 1 tends to zero, and, since H is continuous xP(T > H x ) — ► 0. 



□ 



4 Examples 

Remark 13. The choice F 



and h increasing ca n easily be see n to reduce to the con- 
Valloisl l|l992h with decreasing functions, jvalloisl l|l992i) also has a construction 



struction of 

with increasing functions; this is not a special case of the current construction.) 
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Remark 14. Taking F(s) = h(s) = sj\f2 (or indeed any increasing function of s), we see 
that R(s) — 1 and so H(s) = s. Consequently we have G(l) — I. 1|17|) reduces to 



x 



X t — B t — L t 

which is known to have solution = sup s<t i3 s , and X t is a reflecting Brownian motion 
on (—00, 0]. In particular for some functions <fii,4>2, our stopping rule is now of the form 

T = inf{f > : X t < </>i(L? )} = inf jt > : B t < cj> 2 ^sup B, 

the second form being notable as the same form as the Azema-Yor stopping time; it is easy 
to check the functions are identical. 

We finish by demonstrating how new embeddings can be constructed using the previous 
results, often with useful optimality properties. The embedding we construct will maximise 
the expected value of a convex function of the local time at a level x 7^ 0. (The case where 
x = is the Vallois construction of Remark El) We assume for convenience that x > 0. 
Then we consider the embedding constructed with F(s) — s Ax, and h(s) = s. The resulting 
embedding uses Azema-Yor to embed the left tail of the distribution, until the moment that 
the maximum first reaches x, and then uses the Vallois construction to embed the rest of 
the distribution. The choice of the Azema-Yor embedding will not affect the optimality of 
the construction, but fits in well with the construction described in this paper. We note, 
however, that the construction is the unique construction that maximises E^E'(iy) for a 
convex function and that in this instance L^, = (sup t<T B t ) A x + Lj,. 

Theorem 15. Suppose x > 0. The stopping time T defined in ill'l'l) with F(s) = sAi, and 
h(s) = s maximises E$(Lj.) for any convex function ^> , over all embeddings r of /i in B t 
such that Bt/\ T * s UI. 



Proof. By the optimality of the Vallois construction (see Valloisl |l992lu . the problem 



can 

be reduced to a question of determinining which parts of /i we must embed before we hit 
x. For simplicity, we shall shift the problem to considering the local time at 0, subject to 
having started at some point below 0. This leaves us having to embed the part of /x lying 
above and some mass below dominated by /1. We shall also suppose that /i has a density 
on R given by f(y), however it is clear that this assumption can be removed by considering 
a limiting argument. Consequently, / satisfies 

f(y)dy = l,J yf(y)dy<0. 

We shall prove the theorem using a Lagrangian argument. Suppose the optimal choice 
of embedding has a joint law for stopped distribution of (Bt, Lt) given by p(b, I). Then the 
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problem of finding the optimal embedding can be restated as: 



oo poo 



Maximise / / ^(l)p(b,l) dbdl 

JO J -oo 







subject to / p(b, I) dl = /(&) b>0 (24) 

p(b,l)dl< f(b) b<0 (25) 

p(b,u)dudl (26) 



o 

oo -i poo poo 



J°° bp(b,l)db = ±J_ 



OO J I 



'■ s ^ poo poo 



bp(b,l)db = -- / p{b,u)dudl (27) 
> 2 J_ 00 J l 

P(b,l)>0. 

The necessity and su fficiency of the conditions (|26|) and (|27|l follow from Theorem LI of 



Rovnette et al.l l|2fl02l) . We shall construct the dual problem, and show that, for a suitable 
choice of the Lagrangian multipliers and the conjectured p(b,l), the primal and dual have 
the same value. 

We introduce the multiplier A (b) corresponding to (|24|l for b > and H250 for b < 0, and 
7(Z), rj(l) corresponding to Q26JI and (|27|l respectively. The resulting Lagrangian is: 



p(b,l) 



V(l)-A(b)-i(l)bl {b>0} -r){l)bl {b<0} + l I ( 7 (u) - r)(u)) du 



o 



dbdl 



X(b)f(b)db. 



The corresponding dual problem is therefore to minimise J \(b)f(b) db subject to the term 
in the square bracket being non-positive, and X(b) > for b < (as a consequence of the 
inequality in lj25|0 . For simplicity, we assume also that ^(0) = \&'(0) = 0; the second 
assumption is possible since KLj, = 2K(Bt — x)+, which is fixed for any UI embedding, and 
^ is convex. 

The choice of the Lagrangian multipliers will depend on the optimal construction. As 
such, we are able to conjecture the solution: let 5 < satisfy: bf(b) db = 0. We then 
use the Vallois construction to embed p restricted to [5, oo): let a(l) > decreasing, and 
(3(1) < increasing be the Vallois boundaries, so that T = inf{£ > : B t ^ (f3(L t ), a(L t ))} 
embeds p, restricted to [8, oo), when there is mass f(b) db at 0. 
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Then define 
A(m) 

7(0 



1 1 

+ 



dk; 



l e A(u) f-u 



a ( u ) Jo 
I p A(u) ru 



tf"(m)e~ A(m) dmdu; 



= / ~al — r / *"(™)e (m) dmdu; 
Jo P(u) Jo 



X(b) = 



If 1{b) f u *"(m)e- A Wdmdu : b > 



Note in particular that A(6) > for b < 0. With these definitions, and using the properties 
of the functions a(l) and (3(1), it is now possible to show that 

(0 - \(b) - j(l)bl {b>0} - v(l)bl {h< o } + \J (7(«) - »?(«)) du < 0, 
so that we have a feasible solution to the dual problem, and also 

X(b)f(b) db = m{L T ) 

when T is the candidate embedding. Consequently the embedding is indeed optimal. □ 



/ 
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